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Quantum escapes of two particles with Coulomb interactions from a confined one-dimensional region to a 
semi-infinite lead are discussed by the probability of particles remaining in the confined region, i.e. the survival 
probability, in comparison with one or two free particles. For free-particle systems the survival probability 
decays asymptotically in power as a function of time. On the other hand, for two-particle systems with Coulomb 
interactions it shows an exponential decay in time. A difference of escape behaviors between Bosons and 
Fermions is considered as quantum effects of identical two particles such as the Pauli exclusion principle. The 
exponential decay in the survival probability of interacting two particles is also discussed in a viewpoint of 
quantum chaos based on a distribution of energy level spacings. 
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I. INTRODUCTION 

The escape is a behavior of open systems in which mate- 
rials move out from an observed area. It has drawn consider- 
able attention in various points of view, for example, Kramers' 
escape problem a-decaying nucleus |U@1, the first- 

passage time problem Jll-lstl, the recurrence time problem (01, 
the controlling chaos ff?]. and the Riemann hypothesis js^, etc. 
Escapes involve trans por t, and can be used to calculate trans- 
port coefficients |[l0l - [l2fl . Particles escaping from thermal 
reservoirs can sustain flows such as electric cuiTents iflsifTill . 
Escape phenomena have been investigated in many different 
systems, e.g. billiard systems (by theories |T5l - [l7ll and by ex- 
periments llisi [l9tl ). map systems 1 7 , [sl l20ll . wave dynamics 
II21II22I1 . and stochastic systems ||jl-3], etc. 

A typical quantity to characterize a particle escape is the 
probability of particles remaining within the observed area 
from which particles can move out, the so-called survival 
probability. The survival probability decays in time because 
particles keep to escape from the observed area without com- 
ing back, and its decay properties have been an important sub- 
ject in chaotic dynamics iflll [l2i I23I1 . In classical billiard sys- 
tems, it is conjectured, based on an ergodic argument, that the 
survival probability decays exponentially for chaotic systems, 
while it shows a power decay for non-chaotic systems ifTsll . 
This conjecture has been examined in detail, e.g. in a finite 
size effect of holes lfl6ll . weakness of chaos 101, a connection 
to correlation functions IitIi , and a deviation from an escape 
rate estimated by the natural invariant measure iSHt]. 

Particle escapes have been also discussed in quantum sys- 
tems by using the survival probability. Ref. |24|| discussed an 
escape behavior of a free particle in a one-dimensional system 
by a concrete calculation of wave-packet dynamics, and Ref. 
16(] investigated escapes of a particle with a potential barrier 
Quantum escapes have been also considered by u sing random 
matrix approach for quantum scattering systems l25ll26ll . nu- 
merical approaches to wave-packet dynamics in quantum bil- 
liard systems ll27ll . in a viewpoint of chaotic dynamics. These 
studies show that the survival probability decays in power or 
exponentially, depending on how quantum states are super- 
posed initially, but a general quantum mechanical decay be- 
havior of the survival probability based on chaotic features 



has not been clearly established yet. 

The principal aim of this paper is to investigate many- 
particle effects in escape behaviors of quantum chaotic sys- 
tems in comparison with non-chaotic free-particle systems. To 
investigate them in systems as simple as possible, we consider 
particle escapes from a confined one-dimensional region to a 
semi-infinite one-dimensional lead. Furthermore, as a simple 
chaotic many-particle system we choose the system consist- 
ing of two particles with Coulomb interactions, and discuss 
many-particle effects by comparing the chaotic two-particle 
cases with the non-chaotic cases of one or two free particles. 
In these situations we consider particle escapes whose initial 
states are represented as an energy eigenstate of particles con- 
fined spatially at the initial time, so that voluntariness of initial 
superposition of quantum states in the escape dynamics does 
not appear in discussions of this paper. We calculate survival 
probabilities of such systems as a function of time, and show 
that the survival probability of the Coulomb-interacting two- 
particle system decays exponentially in time, while for the 
free-particle systems it decays in power as a function of time. 
We also discuss quantum effects of identity of two particles, 
like the Pauli exclusion principle, in quantum escape prob- 
lems, appearing as a difference of escape behaviors between 
Bosons and Fermions. It is also shown that a confined system 
consisting of two particles with Coulomb interactions has a 
repulsive feature of energy level spacings, which is regarded 
as a character of quantum chaos. Noting that one-dimensional 
one-particle systems, as well as two free particles in a one- 
dimensional space, cannot be chaotic, this may imply a re- 
lation of the exponential decay of survival probability with 
quantum chaos. 



II. ESCAPE OF MANY PARTICLES IN A SEMI-INFINITE 
ONE-DIMENSIONAL SPACE 

In this paper, we consider quantum systems consisting of N 
particles in a one-dimensional semi-infinite region [0, +cxd). 
Before the initial time f < 0, we set the infinite potential bar- 
rier in the region [I, +00), and confine the particles in the fi- 
nite region [0, /] with a positive constant I. At the initial time 
t we remove this infinite potential barrier in a; > /, so 
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FIG. 1: Schematic illustration of a particle escape in a semi-infinite 
one-dimensional system, (a) Particles confined inside a finite region 
[0, 1] att = 0. (b) Particles escaping to a semi-infinite region at t > 
0, by removing the infinite potential barrier in the region [l,+oo). In 
this situation we call the region [0, 1] as the subspace, and the region 
[l, +oo) as the lead. 



that a particle escape to the region [I, +oo) becomes to occur. 
The schematic illustration of this escaping behavior is shown 
in Fig. [T] (Here, the particles in Fig. [T]are drawn as parti- 
cles with a nonzero finite size to make them visible, but in the 
actual models used in this paper we consider material parti- 
cles with a infinitesimally small size.) To make a clear image 
of this kind of escape phenomena, we call the region [0, 1] as 
the "subspace", and the region {I, +oo) as the "lead," so the 
particle escape occurs from the subspace to the lead. 

This system is described by the wave function 
^{xi, X2, • • • , XNjt) at time i as a solution of the Schrodinger 
equation 



d^{xi,X2 

ih 



dt 



H^{xi,X2, ■ ■ ■ ,XN,t), (1) 



where H is the Hamiltonian operator, h is the Dirac constant 
H, and Xj is the position of the j-th particle, j = 1, 2, • • • ,N. 
If the system consists of identical particles in the quantum me- 
chanical sense, then the wave function 'i{xi,X2, ■ ■ ■ , xat, t) 
must satisfy 

4'(xi,a;2, • • • ,XN,t) 

= ±^'(a;i,X2, • • • ,XN,t)\i^^^^^^^^^^^^^^^ (2) 

for exchange of any particle indices j and k, j = 1, 2, • ■ • , A^, 
k = 1,2, ■ ■ ■ , N . Here, the sign — (+) in ± of the right-hand 
side of Eq. (|2|i should be taken for Fermions (Bosons) ll28ll . 
imposing the Pauli exclusion principle for Fermions. Using 
the wave function '^{xi, x^, - ■ ■ , xn, t) satisfying Eq. ([T]i we 
introduce the probability P{t) defined by 



P[t) = / dxi I dx2 • • ■ / dxN 
Jo Jo Jo 



10 Jo 
x|^'(a;i,a;2,--- ,XN,t)\^- 



(3) 



This is the probability with which A'^ particles are still inside 
the subspace and have not escaped to the lead at time t yet, 
and we call it the "survival probability" hereafter in this paper 
ll29ll . The purpose of this paper is to discuss an escape be- 
havior of particles by a time-dependence of the survival prob- 
ability P{t). It may be noted that the subspace can be re- 
garded as an open system coupled to a semi-infinite lead, but 
the phenomena considered here are not scattering phenomena 
described by a response of the system to incident waves, in 
the sense that particles are always exist only in a finite region 
at any finite time and the wave function is normalizable, i.e. 

J^°° dxi J^°° dx2 - ■ ■ Jq°" dxN |^'(a::i,a;2, • • • ,XN,t)\'^ = 
1 at any time t, different from scattering states including an 
incoming plain wave from the infinite region. 

In general, the survival probability P{t) depends on the ini- 
tial condition. As an initial condition at the time t = 0, in this 
paper we choose an energy eigenstate <i>„(a;i, X2, - ■ ■ , xn) of 
N particles confined inside the subspace region [0, 1] corre- 
sponding to energy eigenvalue En {Ei < E2 < • • • ). In this 
way, we obtain the survival probability for each initial wave 
function ^'(a:;i, X2, ■■• 1 2; jv, 0) = ^n{xi,X2,- ■ ■ ,xn), and 
present it as Pn{t), n = 1,2, ■ ■ ■ . 

Under such initial conditions, we calculate the survival 
probability Pn{t) analytically for free particle cases in Sees. 
Uni and |IV] in this paper. We also show numerical results of 
the survival probability by discretizing space and time in the 
Schrodinger equation for two-particle systems in Sec. |Vl as 
well as for a part of one-particle case in Sec. |III] As a numeri- 
cal technique we use the pseudo-spectral method ll38l[39ll . As 
an example, in Appendix lAl we outline a spatial discretiza- 
tion of the Schrodinger equation for two-particle cases and its 
time-discretization by the pseudo-spectral method, which are 
used to calculate the numerical results shown in this paper 
In these numerical calculations we take the unit of to = 1 
for the particle mass, I = 1 for the length of the subsystem, 
and h = 1 for the Dirac constant. For numerical calcula- 
tions the one-dimensional space is discretized by the length 
Sx = I /JVq with the integer site number Ao of the subspace 
(See Appendix I A II ). For numerical calculations by pseudo- 
spectral method, we also discretize the time by St (See Ap- 
pendix lA2] ). and choose the concrete value of dt so that the 
average energy and the normalization of wave function are al- 
most conserved during the numerical calculations. The total 
system length consisting of the subsystem and the lead in our 
numerical calculations is chosen as L = JVSx = Ml /Mo with 
the total site number A/" (so that the site number of the lead is 
given hy M — Mo), and we calculate the particle escape dy- 
namics in a time interval in which particle's returning back 
to the subsystem from the lead region is negligible. Concrete 
values of the parameters Mo, M, and St will be shown for each 
numerical result in this paper. 



III. ESCAPE OF ONE FREE PARTICLE 

We first consider the case of a single free particle in a 
semi-infinite one-dimensional space, whose Hamiltonian is 
given by _ff = ^[^^/ {2m)]d^ / dx^ with the particle position 
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FIG. 2: Survival probabilities Pn{t),n = 1, 2, 3 of one free particle 
in a semi-infinite one-dimensional space as a function of time t (the 
main figure as log-log plots for a long time behavior and the inset 
as linear-log plots for a short time behavior), corresponding to the 
energies Ei (circle), E2 (triangle), and Es (square), respectively. 
The lines are plots of Eq. (|4} with the coefficient ^ for the cases of 
energy En, n = 1, 2, 3, which are proportional to t^"^. 



X = xi. Using this Hamiltonian we solve the Schrodinger 
equation ihd'^ {x ,t) / dt = H'^{x,t) for the wave function 
^'(a;, t) of the system for a; > and ^"(0, t) = 0, and calcu- 
late the survival probability (|3]l. 

In this system, as shown in Appendix IB II we can solve the 
Schrodinger equation analytically, and the survival probability 
P{t) is represented asymptotically as 



Here, the constant Ai is given by 



2 / ml 



dx X '^{x,{)) 



(4) 



(5) 



with the wave function '^{x, 0) at the initial time t = 0. 
Equation (|4]l means that the survival probability P{i) decays 
by power ^ t^^ asymptotically in time, for arbitrary initial 
conditions of the subspace as far as Ai ^ 0. Power decays 
of the survival probability for one-particle systems in a one- 
dimensional space have been discussed in some papers ll6l l24ll . 

Figure |2] is the graphs of the survival probabilities Pn{t), 
n = 1,2,3 obtained by solving the Schrodinger equation for 
the one free particle numerically, using the subspace site num- 
ber TVq = 60, the total space site number TV = 32768, and 
the discretized time interval 5t = 10^^. The energy values 
corresponding to these graphs are i?i = 4.77, i?2 = 19.1, 
and ~ 42.9. For a comparison, in this figure we also 
draw the graph of Eq. (|4]i with the coefficient ^ for each 
energy. The numerical results of the survival probabilities 
Pn{t),n = 1,2,3 in Fig. |2]show a behavior of the power 
decay ~ t^^ independent of n in a large time region, and they 
are consistent with our analytical result Q including value of 
the coefficient (|5]l. 



FIG. 3: Probability distribution function f:i{x,t) = \^3{x,t)\^ of 
the particle position as a function of t and x, corresponding to the 
energy iJg, for one free particle in a semi-infinite one-dimensional 
space X > 0. 



In order to visualize an escape behavior of one-free-particle 
systems, we show in Fig. [3] the probability distribution func- 
tion fn{x,t) = |'I'„(a;,<)p for ?i = 3 as a function of time 
t and position x. Here, we took the case of n = 3 in Fig. [3] 
so that this spatial distribution fsix, t) has three peaks at the 
initial time t — 0, but two of these three peaks decay quickly 
in time and only one peak survives for a long time and moves 
away from the subspace region [0,1]. 



IV. ESCAPE OF IDENTICAL TWO FREE PARTICLES 

As a many -particle effect in escape phenomena, we first dis- 
cuss the quantum effect of identity of two particles, such as the 
Pauli exclusion principle, in a semi-infinite one-dimensional 
system. The Hamiltonian operator of the system is given by 
H ^ -[ftV(2m)] (d^/dxl + d^/dxl) with the position Xj 
of the j-th particle {j = 1, 2), and we impose the condition 
^ for the wave function X2, i) as 



^{xi,X2,t) = ±'il{x2,Xi,t). 



(6) 



Here, the sign + (— ) in the right-hand side of Eq. (|6]l is taken 
when the particles are identical Bosons (Fermions). The con- 
dition (|6]l is guaranteed at any time t as far as it is imposed at 
the initial time t — because the Hamiltonian operator H is 
symmetric for exchange of two positions of the particles. 



A. Boson Case 

In the case of identical two free Bosons, as shown in Ap- 
pendix IB 2 al the asymptotic behavior of the survival proba- 
bility P{t) is represented as 



Pit) 



A 



2b 



(7) 
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Here, is defined by 
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/ dxi / dx2 XiX2^{xi,X2,0) 

Jo Jo 



(8) 



with the initial wave function ^(xi, X2, 0). Therefore, in the 
case of identical two free Bosons, the survival probability 
P{t) asymptotically decays in power ^ t^^, which is sim- 
ply the square of the result of the one free particle discussed 
in Sec. Hn] 




B. Fermion Case 

In the case of identical two free Fermions, as shown in Ap- 
pendix IB2b[ the asymptotic behavior of the survival proba- 
bility P{t) is represented as 



Pit) '^^^ ^ 



(9) 



with the constant A2 / defined by 

10 

'/ / m l. \ 

A2f 



2 /ml 
47257r2 It 



dxi / dx2 



XX1X2 {xl - xl) 'i'{xi,X2,0) 



(10) 



It is important to note that in the case of identical two free 
Fermions, the survival probability P{t) asymptotically decays 
in the power ^ which is different from the power ^ t^^ 
in the corresponding Boson case shown in Eq. (|7]i. In the 
identical two free Fermions the power decay term ~ t of the 
survival probability disappears by an cancellation for the Pauli 
exclusion principle ^{xi,X2,t) + ^{x2,xi,t) = 0, and such 
Fermions escape faster qualitatively than the corresponding 
identical two free Bosons. 



V. ESCAPE OF TWO PARTICLES WITH COULOMB 
INTERACTIONS 

Now, we discuss an escape behavior of a semi-infinite-one- 
dimensional system consisting of two particles with Coulomb 
interactions. The Hamiltonian operator of the system is given 

hyH = -[^7(2™)] {d^ / dxj+d^/ dxi) +Uixi,x2) where 
U{xi , 2:2) is the Coulomb potential energy as 



U{xi,X2) 



A 



v/rf2 + (xi - X2y 



(11) 



Here, A is given by A = g^/(47reo) with the particle charge 
q and the dielectric constant eq, and d is a small but nonzero 
constant appearing as an effect of a quasi-one-dimensionality 
of the system OOll . Using this Hamiltonian we solve the 
Schrodinger equation ift9\E'(a;i, X2, = H'i/{xi,X2,t) 



FIG. 4: Survival probabilities P„{t),n — 1, 2, 3 for two particles 
with Coulomb interactions in a semi-infinite one-dimensional space 
as a function of time t (linear-log plots), corresponding to the ener- 
gies El (circle), E2 (triangle), and E3 (square), respectively. The 
inset is linear-log plots of Pn{t) for these energies in a short time 
region. The lines are the fits for each survival probability P„ (t) to an 
exponential function a exp(— with fitting parameters a and /3. 



for the wave function '^{xi,X2, t) of this system with the con- 
dition (|6]l, then calculate the survival probability (O. 

In this section, we show only results of the identical-two- 
Fermion case, since the corresponding results for the Boson 
case are quite similar to the Fermion case and we could not 
find a particular difference (e.g. in an exponential decay be- 
havior of the survival probability and value of its escape rate 
as will be shown for Fermions in this section) between the 
Fermion cases and the corresponding Boson cases in our nu- 
merical results. 

Figure|4]is the survival probabilities Pn{t), n = 1, 2, 3 at 
time t for identical two Fermions with Coulomb interactions 
in a semi-infinite one-dimensional space. For numerical cal- 
culations to obtain the graphs in Fig. HI we used the parameter 
values A = 60 and d = 10^''/6 for the potential energy ( fTTT ). 
and also the subspace site number A/q = 60, the total one- 
dimensional space site number M ~ 2048 (so the total site 
number in two-dimensional X1X2 space is M"^ = 4194304), 
and the discretized time interval 5t = 10^"^. The energy 
values corresponding to the survival probabilities P„ {t),n = 
1,2,3 in Fig. Hare given by Ei = 149, E2 = 202, and 
i?3 = 213, respectively. By the repulsive Coulomb interac- 
tion, the survival probability P„ {t) decays much more rapidly 
than the case of one free particle. Figure |4] shows that the de- 
cay of the survival probability P„(i) for Coulomb-interacting 
two particles is well approximated as an exponential decay 
after a short time, different from the cases of one particle dis- 
cussed in Sec. |III] In Fig. |4]we also showed fits for each sur- 
vival probability f „ {t) to an exponential function a cxp(— /3t) 
with the fitting parameter a and p. Here, the values of fit- 
ting parameters are chosen as (a, /3) = (5.03 x 10^, 166) for 
Pi(<), {a, 13) ^ (4.81 X 10^167) for P2{t), and [a, (3) = 
(3.21 X 10^, 176) for P:i{t). The escape rate, defined as the 
parameter /3, does not depend strongly on for value of the en- 
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FTG. 5: One-particle spatial distribution function fs{x,t) = 
ff^°° dx2 1^3(2;, X2,t)\'^ as a function of time t and position x, cor- 
responding to the energy £3, for two particles with Coulomb inter- 
actions in a semi-infinite one-dimensional space a; > 0. Two peaks 
are seen in the time region t > 0.1, corresponding to two particles. 




FIG. 6: Distribution function D[s) of energy level spacings (solid 
line) of two particles confined in the one-dimensional region [0, 1] 
without lead. The broken line and the dotted line correspond to the 
Wigner distribution and the Poisson distribution, respectively. 



ergy En in our numerical results. 

Figure |5] is a graph of the one-particle spatial distribu- 
tion function fn{x,t) = J^°° dx2 |5'„(a;, X2, t)P = 

J^°° dxi |^'„(a;i, x, t)p for n = 3 as a function of time t 
and position x. In this figure it is clearly seen that two peaks 
corresponding to two particles can survive for a long time al- 
though there are many peaks in a short time region. The posi- 
tion of one peak of these two peaks closer to the wall at a; = 
does not move very much in time, probably because the par- 
ticle corresponding to this peak receives a repulsive Coulomb 
force from another particle to the direction of the wall and its 
escape movement is suppressed. On the other hand, the posi- 
tion of another peak moves away quickly from the subspace. 

VI. CHAOS IN IDENTICAL-TWO-PARTICLE SYSTEMS 
WITH COULOMB INTERACTIONS 

Now, we discuss a difference between one free particle and 
Coulomb-interacting two particles in a one-dimensional space 
in a viewpoint of quantum chaotic dynamics. 

Chaos is defined dynamically as a classical system with 
a dynamical instability, i.e. a strong sensitivity to a small 
difference of initial conditions leading to at least one posi- 
tive Lyapunov exponent, and quantum chaos has been inter- 
preted as a quantum mechanical system whose classical coun- 
terpart is chaotic ll3ll - [33ll . In this sense, there is no quantum 
chaos in any one-dimensional one-particle system with any 
time-independent potential, because the corresponding clas- 
sical system does not have the dynamical instability. On the 
other hand, quantum two-particle systems in one-dimensional 
space can be chaotic, because the coiTesponding classical sys- 
tem can have a dynamical instability. The quantum chaos has 
been also widely investigated by the distribution of energy 
level spacings ll3ll - t33ll . Following to this idea, the system with 
the time-reversible dynamics is regarded chaotic in a quan- 



tum sense if the distribution of energy level spacings is close 
to the Wigner distribution (7r/2)x exp(— 7r.T^/4) by the Gaus- 
sian orthogonal ensemble (GOE), while it is non-chaotic if the 
distribution of energy level spacings is close to the Poisson 
distribution exp(— a;). Energy level spacing distributions of 
many-particle systems with Coulomb interactions have been 
calculated for two-particle cases in a one-dimensional space 
llsoll . in a one-dimensional space with a random potential i34ll . 
and in a two-dimensional space |[35i[36ll . as well as for three- 
particle cases IstIi . 

To calculate the energy level spacings of a quantum sys- 
tem we have to eliminate some trivial symmetric degrees of 
freedom from the Hamiltonian of the system. The Hamilto- 
nian of identical-two-particle systems with the Coulomb po- 
tential ( fTTl i is invariant for exchanging the two particle posi- 
tions xi and X2- To eliminate the degeneracy caused by this 
symmetry from energy level spacings we calculated the en- 
ergy eigenvalues of the spatially discretized Hamiltonian ma- 
trix (shown in Appendix I A 11) only for xi > X2- (Here, we 
used the Fermion case so that by the Pauli exclusion princi- 
ple the possibility of xi = X2 for the two particle positions 
xi and X2 is zero.) Further we also separate this Hamilto- 
nian operator into the symmetric part and the anti-symmetric 
part for the transformation xj — > / — Xj,j = 1,2, and we 
take into account only of its symmetric part. To calculate the 
distribution function of energy level spacings, we first calcu- 
late the energy eigenstates £„, n = 1, 2, • • • , /i, £1 < £2 < 
• • ■ < S/j. of this reduced Hamiltonian. Then, we fit the in- 
dices n — /io, /^o + 1, • • • , /i as a function of energy £„ by a 
smooth function p{£) = X]fe=i o-k£'' with fitting parameters 
flfe, fc = 1, 2, • • • , 10. Here, we did not take into account the 
energy with small indices n = 1, ■ • ■ , ^0 ^ 1 for this cal- 
culation with an integer yUo > 1- in order to exclude strongly 
quantum regime from calculation of energy level spacings. 
Using this fitting function we calculate the distribution func- 
tion D{s) of energy level spacings as the distribution function 
of s = p{£j+i) - p{£j), j = po,po + l, - ■ ■ ,p-l. 
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Figure |6]is a graph of the distribution D{s) of energy level 
spacings (solid line) for a Coulomb-interacting identical-two- 
particle system confined in a one-dimensional region [0, /] 
without lead. Here, we used the same values of the system 
parameters m, h, I, A, d, and TVo as used in Sec. [V] and also 
/io = 20. For comparisons, in Fig. |6]we also show the Wigner 
distribution (broken line) and the Poisson distribution (dotted 
line). The distribution D{s) in Fig. |6]shows a repulsive behav- 
ior of energy level spacings and is different from the Poisson 
distribution, although its repulsion is not strong enough for the 
distribution D{s) to be fit to the Wigner distribution. In this 
sense, the identical-two-particle system with Coulomb inter- 
actions in a confined one-dimensional space is regarded to be 
weakly chaotic. 



VII. CONCLUSION AND REMARKS 

In this paper we discussed escape behaviors of one- and 
two-particle systems from the one-dimensional finite region 
(0, /) to the semi-infinite one-dimensional lead [I, +oo) with a 
positive constant /. We prepared the initial condition of wave 
function of the system so that at the initial time t = the 
wave function is given by an energy eigenstate of the particle 
system confined in the finite region (0,1). Under these ini- 
tial conditions we calculated the survival probability for par- 
ticles to stay within this region (0,1) at time t. We showed 
that the survival probability decays in power ~ tr^ asymp- 
totically in one-dimensional one-free-particle systems, and it 
decays exponentially in two-particle systems with Coulomb 
interactions. The quantum effect of identity of many particles 
such as the Pauli exclusion principle in a behavior of the sur- 
vival probability were also discussed in the cases of Bosons 
and Fermions, and it was shown that the survival probability 
decays asymptotically in power ~ for the identical two 
free Bosons but in power ^ t~^^ for the identical two free 
Fermions. 

This study is motivated to investigate, not only quantum and 
many-particle effects, but also a chaotic effect in escape phe- 
nomena. In classical mechanical cases there is a conjecture 
by which the survival probability decays exponentially for 
chaotic systems while it decays in power for non-chaotic sys- 
tems. To check a chaotic feature of two-particle systems with 
Coulomb interactions we calculated the energy level spac- 
ing of the system in a confined one-dimensional space, and 
showed a repulsive feature of energy level spacings in this 
system. This result suggests that two-particle systems with 
Coulomb interactions in a confined one-dimensional space are 
weekly chaotic. On the other hand, one-free-particle systems 
in a one-dimensional space is not chaotic, because its corre- 
sponding classical systems do not have a dynamical instabil- 
ity. Therefore, our results for quantum escapes in this paper 
may be consistent with the classical conjecture for a chaotic 
effect in exponential decays of survival probabilities. 

As a remark, although two-particle systems with Coulomb 
interactions show an exponential decay in the survival proba- 
bility and they are (weekly) chaotic, it might not be necessary 
to mean that chaos is the origin of this exponential decay. In 




FIG. 7: Survival probabilities P„{t),n = 1,2,3 of a one-particle 
system with a localized single impurity potential in a semi-infinite 
one-dimensional space as a function of time t (linear-log plots), cor- 
responding to the particle energies Ei, E2, and E3, respectively. 
Here, the impurity is located just outside the subspace, and mag- 
nitude of the impurity potential is chosen to be very strong in com- 
parison with the particle energies. 



other words, we still should be careful to a possibility that a 
chaotic feature is not a necessary condition for an exponential 
decay of survival probability in quantum systems. As a result 
related to this point, in Fig. |7] we show the survival proba- 
bilities Pn{t),n = 1, 2, 3 for the one-particle system with a 
very localized single impurity at the just outside of subspace 
in a semi-infinite one-dimensional space, where we used pa- 
rameter values as the subspace site number A/q = 60, the total 
space site number A/" — 16384, and the discretized time inter- 
val St = 10~^. Here, we chose the magnitude of the impurity 
potential as 1.8 x 10'^, which is much larger than the parti- 
cle energies Ei ~ 4.77 (circle), E2 ~ 19.1 (triangle), and 
Es ~ 42.9 (square) in Fig. [7] The one-particle system in a 
one-dimensional space is not chaotic, but Fig. |2]shows a clear 
exponential decay in the survival probability of such a one- 
particle system with a strong impurity potential. However, in 
Fig. |2]the exponential decay rates of the survival probabilities 
for this one-particle case with an impurity strongly depend on 
the value of their particle energies, different from results of 
the two-particle systems with Coulomb-interactions as shown 
in Fig. m 

It may also be meaningful to mention that an exponential 
decay of the survival probability in interacting two-particle 
cases shown in this paper occurs in a rather strong Coulomb 
interaction regime. In cases with much smaller interaction 
magnitude A than used in Sec. |V] we observe a decay close 
to be in power, which is similar to the asymptotic decay be- 
havior for identical two free particles shown in Sec. |IV] In 
a system with such a small A, the distribution of energy level 
spacings is rather close to the Poisson distribution, implying 
that the system is rather close to be non-chaotic. This may be 
a supportive evidence that a chaotic feature may be related to 
an exponential decay behavior of the survival probability at 
least in two-particle systems. 



7 



One may notice that in this paper exponential decays in 
two-particle systems with Coulomb interactions were shown 
numerically, so there is still a possibility that it might be only 
be a finite time property and it might decay in power in the 
long time limit t — > +00. Related to this remark, Ref. |0] ar- 
gued that the survival probability for one-particle systems in 
an effectively one-dimensional space with a delta-functional 
impurity shows a power decay, even if it decays exponen- 
tially in time for a finite time region like in Fig. [7] How- 
ever, it would be valuable to show some results about dif- 
ferences between one-particle (non-chaotic) cases and two- 
particle (chaotic) cases even if they are finite time properties. 

As a future problem on the subject of this paper it would 
be interesting to investigate escape behaviors of systems con- 
sisting of more than two particles, although in order to study 
such a large system numerically we would need much better 
numerical resources and techniques than used to obtain re- 
sults in this paper Such large systems with interactions could 
be more strongly chaotic than two-particle systems, whose 
chaotic strength was not so strong as shown in Sec. [Vll and 
they could give better situations to study quantum chaotic ef- 
fects in escape phenomena. We could also consider a par- 
ticle escape of such a large number of particle systems as 
a driving source to produce a particle current from a parti- 
cle reservoir. For such a case, as the initial state we could 
choose an equilibrium (e.g. canonical) state, which is repre- 
sented as the density matrix produced by the energy eigen- 
states {$„(a;i, X2, • ■ • ,XN)}n used in this paper with the 
weight of an equilibrium distribution. 
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Appendix A: Numerical Calculations of a One-Dimensional 
Two-Particle Dynamics 



semi-infinite one-dimensional space by a positive constant SI, 
= 0,1,2,- •• for j = 1,2. In this dis- 



SO Xi 



rijSl, Tij 



cretization of the one-dimensional space, the spatially second 
derivative d^/dx^ applying to any function X{xj) is repre- 
sented as 

d'^X{xj) X{nj + 1) - 2X{nj) + X{n.j - 1) 



dx] 



(Al) 



k=0 



where X{nj) = X(nj5l) is the spatially discretized function 
of X{xj). Here, the matrix /C = (ICjk), j = 0,1,2,---, 
k ~ 0,1,2, - -, is the matrix whose only nonzero elements 
are % = 2 and /Cq+i)j = = -1- J = 0, 1, 2, • • • . 

Using Eq. dAlb . the Hamiltonian operator ti is represented as 
the matrix H for the spatially discretized representation: 



H ^ H = 



2m5V 



;K + U{n)I 



(A2) 



with n = (711,712) and f/(n) = U{niSl,n2Sl) using the po- 
tential U{xi^X2) given by Eq. (fTTT t for the continuous space 
case. Here, the matrices K = (A'nnO and / = (/nn') are 
defined by 



nn' — n'^ ^n^n^ 

/nn' 



(A3) 
(A4) 



respectively, for any n = (ni,n2) and n' = {n'^,n'2)- Using 
Eq. ( IA2b the Schrodinger equation is spatially discretized as 
ihd'^{t)/dt = H'^{t) as a equation for the vector = 
($(n, i)) defined by l>(n, t) = '^{niSl, n25l, t) with the vec- 
tor index n. 

It may be meaningful to represent the Hamiltonian matrix 
H by using the Dirac notation. Introducing the state |n) as 
the one forming a complete {J2n \^) ~ 1) ^^'^ orthogonal 
((n|n') = Sn^n[Sn^n'^) Set for the site indexes n and n', the 
Hamiltonian matrix H can be represented as the operator 



(A5) 



n.n 

(|n-n'| = l) 



In this Appendix we represent how we discretized the 
Schrodinger equation for the system consisting of two par- 
ticles in a semi-infinite one-dimensional space to calculate the 
survival probability P{t) in this paper Especially, we show 
a spatially discretized Hamiltonian and outline the pseudo- 
spectral method to solve the Schrodinger equation discretized 
in time. 



where £„ and u are defined by e„ = [2h'^ / (mSP)] + 
U{niSl,n26l) and u = —h'^/{2m5P), respectively. The 
operator iA5\ has the same type of form as a tight- 
binding Hamiltonian with the site energy £„ and the hop- 
ping rate u, and the Schrodinger equation is represented 
as ihd\'^{t)) /dt ~ H\'i'{t)) as a equation for the state 
l*W)^E„'5(n,t)|n). 



1. Spatial Discretization of the Hamiltonian Operator 

We consider two particles in a semi-infinite one- 
dimensional space and take Xj{> 0) as the position coordi- 
nate of the j-th particle, j = 1,2. Then, we discretize the 



2. Time-Discretization of Schrodinger Equation by 
Pseudo-Spectral Method 

In the previous subsection of this Appendix we discussed 
how we spatially discretized the Hamiltonian operator. In this 
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subsection we outline how we discretize the time-evolution 
by the Schrodinger equation in the way called by the pseudo- 
spectral method. 

We consider a one-dimensional space of the length L con- 
sisting of the subspace and the lead, and note that the func- 
tion A" (.T 1 , 2:2 ) of .T 1 and X2 satisfying the boundary condition 
X{Q,X2) = X{xi,Q) ^ X{L,X2) = X{xx,L) = can be 
Fourier-transformed as 



^"(^1,^2) = VT / / dx2X{xi,X2) 











X sm| I smi ~Y^^2 



= P[X{xi,X2)] , 



X{xi,X2) 



2 +00 +00 



fel=l fe2=l 



x(k^M) 



'Trfci \ . / 7rfc2 
X sm| I smi -^.T2 



= f-^ X{h,k2) 



(A6) 



(A7) 



by using the relation Jq dx sm{TTkx/L)sm{'Kk'x/L) = 
L5kk'/X k = 1, 2, • • • , fc' = 1, 2, • • • etc. Using Eq. ( lATl i 
we obtain 



KX{xi,X2)=J'^^ K{ki,k2)X{ki,k2) 



(A8) 



with the kinetic operator A' = —[l/(2m)](9^/9xf + 9^/9x3) 
and the function K{ki, k2) = [k'^ / {2mL'^)] {kf + of ki 
and k2 . 

We discretize the time by a positive constant 6t, so t — > 
lySt, 1/ = 0,1,2, ■■■ . By using the formal solutions of the 
Schrodinger equation, the wave function 'i/{xi,X2,t + 6t) at 
time t + 6t.is related to the wave function 'i>{xi,X2,t) at time 
t as 



^{x,,X2,t + St) = e~'™/'"'*(xi,X2,t) 



(A9) 



,-iU(xi,X2)St/(2h) p- 



-iStK(kiM)/fi j: 



e-*^(^i^^^)**/(2ft)*(xi,a;2,t) + 0{6f) (AlO) 



where we used the relation ( IA8I ) and the boundary conditions 

4'(0,X2,t) = *(.Tl,0,t) = *(i,X2,t) = = 

for the wave function of the system used in this paper. By Eq. 
( I A lot we can calculate the wave function '^{xi,X2,t + 5t) at 
time t + 5t from the wave function '^{xi,X2,t) al the previous 
time t. 



An advantage of the pseudo-spectral method dAlOl l is that 
in this method we do not have to apply the space-differential 
operator Q-^HSt/h -^^ j-j^g time-evolution of wave function, and 
it is replaced by simple multiplications of just the numbers 
^-iU(xuX2)St/(2h) jjjjj ^^tStk{kiM)/h^ leading to less com- 
plicate numerical calculations than to use Eq. iA9i directly. 
Instead, we need to do a Fourier transformation and an inverse 
Fourier transformation for one step of the time evolution, but 
we can use the technique called by the fast Fourier transfor- 
mation lf39ll in actual numerical calculations. The fast Fourier 
transformation requires the calculation time proportional to 
be A/'log2 A/" (instead of N"^) for the total {xiX2 space) site 
number A/", and it is a big advantage for a fast numerical cal- 
culation of large AA systems such as used in this paper. 



Appendix B: Survival Probability of Free Particle Systems in 
the semi-infinite one-dimensional space 

In this appendix we calculate the survival probability an- 
alytically for free particle systems in a semi-infinite one- 
dimensional space. First, we calculate it for the case of one 
free particle, and show analytically an asymptotic power de- 
cay (HJi of the survival probability P{t). Secondly, we show 
the different power decay behaviors O and (|9]l of the survival 
probabilities between identical two free Bosons and Fermions 
with the quantum effect of identity of two particles, such as 
the Pauli exclusion principle for Fermions, in a semi-infinite 
one-dimensional space. 



1. One Free Particle Case 

For one free particle in a one-dimensional space, the Hamil- 
tonian operator is given by H = — [fi,^/(2m)]9^/9a;^. Then, 
the wave function ^ {x, t) for this Hamiltonian system in 
the full one-dimensional infinite region (— cxo, +00) is repre- 
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sented as ||4C 



27riht 



dy ^ (y,0) 

im{x — y)^ 



X exp 



2ht 



(Bl) 



for any initial wave function 5* (x, 0) of the system. Using 
this function '5 {x, t), the wave function '^{x, t) for the sys- 
tem with the same Hamiltonian but in the semi-infinite one- 
dimensional region [0, +oo) is given by 



E'it)-'{^'{x,t)-^'i-x,t)] (B2) 

-t-oo 

dy ■^{y,0) 



m 
2-Kiht 



X exp 



i'm(x — yY 
2ht 



(B3) 



for a; > 0, so that the boundary condition ^(0, t) = in the 
hard-wall at x = is automatically satisfied at any time t by 
Eq. dull. Here, = dx {x,t) ~ ^' {-x,t)\^ 

is the quantity to normalize the wave function '^{x^t) as 
^ dx \^[x,t)\^ = 1 for the semi-infinite space [0, +cx3), 
and the initial wave function ^'(a;, 0) = S (0)^^[^' (.t, 0) — 
^'(-x, 0)] satisfies the condition '^{-x,Q) = 0) for 

any real number x in Eq. ( IB3I ). By using Eq. ( IB3I ) and noting 
the fact that value of the initial wave function ^'(a:;, 0) is zero 
for |a;| > I, the wave function ^'(a;, t) of one free particle in 
the semi-infinite region [0, +oo) is represented as 



■^{x,t)^ / dyG{x,y-t)^!{y,Q) 
Jo 

at the position x at time t with G{x, y; t) defined by 

/ Tfl 



(B4) 





im{x — yY 


jexp 


2ht 



exp 



im{x + yY 

m 



(B5) 



as the time-evolutional propagator for one free particle in the 
semi-infinite one-dimensional space [0, +00). 

For large t we expand the propagator G{x, y; t) as 

f2i /"m\3/2 



(i) - + <«« 



G{x,y]t) 



up to the smallest non-zero order of 1 /t. Inserting Eq. 
into Eq. (IB4b and then calculating the survival probability Q 
for = 1, we obtain Eq. (|4]i with the coefficient 



2. Identical Two Free Particle Cases 

We consider an identical-two-free-particle system with no 
potential energy in a semi-infinite one-dimensional space. In 



this case, because of no potential energy, the time-evolution 
of the wave function 5'(.ti, 3:2, t) of this system is dominated 
by the one-particle propagator ( IB5I ) and is given by 



(B7) 



*(xi,a;2,t) = dyi dy2 G{xi,yi;t)G{x2,y2;t) 
Jo Jo 

x*(2/i,y2,0). 



Here, we used the fact that values of the wave function 
^{yi,y2, 0) at the initial time f = are nonzero only in the 
region satisfying < yi < I and < ?;2 < ^ as assumed in 
this paper. 

Now, we impose the condition 



*(a;2,a;i,0) = ±«'(a:i, X2, 0), 



(B8) 



i.e. the condition (|6]l at the initial time t = 0. Here, the sign 
+ (— ) in the right-hand side of Eq. ( IB8I ) is taken for Bosons 
(Fermions). Under the condition ( IB8I ) the wave function ( |B7| ) 
automatically satisfies the condition (|6]l at any time t. We can 
rewrite Eq. ( IB7b using Eq. ( IB8I 1 as 



*(a;i,a;2,t) = - [^ixi,X2,t) ±'^{x2,xi,t)] 



o , dyi dy2 [G{xi,yi]t)G{x2,y2\t) 
^ Jo Jo 

±G(x2, 2/1; t)G(xi, 2/2; t)] 2/2, 0). (B9) 

so that the wave function ^(a;i, X2,t) automatically satisfies 
the condition ^ at any time t without noting the condition 
anymore. 



a. Boson Case 

For the system consisting of identical two free Bosons, 
by Eq. ( IB6I 1 the quantity G{xi,yi-t)G{x2,y2',t) + 
G(x2, yi; t)G(xi, 7/2; t) is asymptotically represented as 



G{xi,yi;t)G{x2,y2;t) + G{x2,yi]t)G{xi,y2]t) 

(BIO) 



4z / TTi \ 



Inserting Eq. (IB 10b into Eq. ( |B9b and then calculating the 
survival probability ([3j for iV = 2, we obtain Eq. Q with the 
coefficient ([8]). 



b. Fermion Case 

For the system consisting of identical two free Fermions we 
expand the one-particle propagator ( IB5b as 



G{x,y;t) 

l2i /m 



3/2 



1 / HIj \ ^ 

--(-) (3.^ + 10a.V + 3y^) 



or 



(BID 
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up to the order of X/t^/"^, which is higher order than in Eq. 
( |E6] l. By Eq. dBllb the quantity G{xt .VT-.t^Glxo.v-y-.t) - 
G{x2, yi\ t)G(xi, 2/2; t) is asymptotically represented as 

G{xi,yi]t)G{x2, y2]t) - G{x2,yi;t)G{xi,y2;t) 
= ^ (^) xiX2yiy2 {xl - xl) {yl - yl) 
+0{t-^) . (B12) 



Inserting Eq. (IB12l i into Eq. ( |B9l l and then calculating the 
survival probability (O for iV = 2, we obtain Eq. (|9]l with the 
coefficient ( fTot . 
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